Aspect-Ratio Scaling and The Stiffness Exponent 9 for Ising Spin Glasses 
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We introduce the technique of aspect-ratio scaling to study the scale-dependence of interfacial en- 
ergies in Ising spin glasses, and we show how one can use it to determine the stiffness exponent 9 in 
a clean way, with results that are independent of the domain-wall-forcing boundary conditions im- 
posed on the system. In space dimension d = 2 we obtain 6 = —0.282(3) for a Gaussian distribution 
of exchange interactions. 



The determination of stiffness exponents in spin glasses 
has a long history, going back to the early 1980s |]-|. 
Loosely speaking, the stiffness exponent 9 of an Ising 
spin-glass is defined by the statement that the energy, 
Eint, of an interface between ground states scales with 
length scale L as Eint ~ . The looseness in this defi- 
nition comes from the vagueness surrounding the phrase 
"length scale L" . Since such interfaces are fractals |,|, 
the length L does not refer to any measure on the inter- 
face itself, but rather to the size of the region in which 
the interface is confined. Traditionally, square or (hy- 
per)cubic regions of side L are employed, with an inter- 
face imposed by a suitable change of boundary condition. 
The problem with this approach is that the results often 
seem to depend on the choice of boundary conditions. 
Here we explain why this is so and introduce the tech- 
nique of aspect-ratio scaling as a method of obtaining 
a well-defined value for 9, independent of the boundary 
conditions imposed. 

First we briefly review some of the different boundary 
conditions which have been used. For simplicity we dis- 
cuss only space dimension d = 2, but generalization to 
d > 2 is trivial. In each case the boundary condition in 
the y-direction is periodic, while the domain-wall-forcing 
boundary conditions are imposed in the a:-direction. 

(i) Periodic- Antiperiodic (P-AP): Ground state en- 
ergies, Ep and Eap, are determined for periodic and 
antiperiodic boundary conditions respectively (the lat- 
ter often defined by reversing the signs of one column 
of bonds). The lower-energy state is obtained for P or 
AP boundary conditions with equal probability, and the 
higher-energy state contains a domain wall relative to 
the lower. The interface energy is therefore defined as 
Eint = \Ep — Eap\- 

(ii) Free-Antifree (F-AF): The ground state energy, 
Ep, and spin configuration are obtained with free bound- 
aries in the x-direction. The spins at one end are held 
fixed, those at the other end flipped, and the new ground 
state energy, Eap, obtained. In this case Eint = Eap — 
Ep, and E,nt > 0. 

(iii) Random- Antirandom (R-AR): The spins at both 
ends are clamped in random configurations, and the 
ground-state energy, Ep, is obtained. The spins at 
one end are held fixed, those at the other end flipped, 
and the new ground-state energy, Ear, is found. Now 
Eint — \Er — Ear\- 



We now discuss the application of these methods to 
the nearest-neighbor Ising spin-glass model with Hamil- 
tonian H = — J2{ij) JijSiSj, restricting our attention ini- 
tially to a Gaussian distribution of the exchange interac- 
tions, Jij, in dimension d ~ 2. The P-AP method has 
traditionally been the most popular. The exponent 9 has 
been measured for d = 2 |§, 3 |,§ and 4 g. For d = 2 
the result 9 = —0.281(2) was obtained from square sys- 
tems of size L < 30 . A recent result by Hartmann and 
Young (HY) 1^ on much larger square systems {L < 480), 
but using free boundaries in the y-direction, is consistent 
with this: 9 = -0.282(2). 

R-AR boundary conditions were used in early work 
by two of us |] to obtain 9 = -0.291(2) for L < 12. 
We believe (as will be discussed further below) that the 
small discrepancy with references is due to the small 
range of sizes available in the earlier study. 

Results for F-AF boundary conditions differ from the 
P-AP results by a somewhat larger amount, with 9 — 
-0.20 found for i < 24 0, and 61 = -0.266(2) obtained 
by HY using L < 320 (but with free boundaries in the 
y-direction). HY conjecture that the exponent is actually 
independent of the boundary conditions, but that larger 
sizes {L ;» 320!) are needed to reach the asymptotic 
regime for F-AF boundary conditions. 

The questions raised by these results are (i) Are the 
results really boundary-condition-independent? If not, 
what does 9 mean, and if 9 is not well-defined, what are 
we to make of the conventional result ^ ~ T^/^ ^ for 
the divergence of the correlation length as T — > 0? (ii) 
If 9 does have a well-defined meaning, independent of 
the boundary conditions, is there an efficient method to 
obtain its asymptotic, boundary-condition-independent 
value? 

To answer these questions we introduce here the idea 
of aspect-ratio scaling (ARS). Using this approach we ob- 
tain strong evidence for a unique 9, and we determine its 
value as 6* = —0.282(3), consistent with values quoted 
above from studies using P-AP boundary conditions. 
For a given a number of spins, this approach apparently 
converges much more rapidly than using square samples. 

Consider a system of length L and width M , where 
we will usually take L > M . The very natural assump- 
tion underlying ARS is that the mean interfacial energy 
(averaged over samples) has the asymptotic form (for L 
and M both large) 
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where L/M = R is the aspect ratio of the samples. Now 
consider the hmit i? — > cxd. In this hmit the system be- 
haves hke a d = 1 system, for which one can show that 
P,p^ (Eint) ~ Imposing this hmiting form on (|l|) 

requires F{x) ^ 1/x for x oo, and gives 



{Eint) 



L 



L:$> M . 



(2) 



It is also of interest to consider the limit M L. In 
this limit, sections of the interface whose spatial extent 
is much larger than M are essentially independent, so 
we expect (for an M"*"^ x L system) the Af-dependence 
(Eint) ~ M^'^-i^/^ for P-AP or R-AR boundary condi- 
tions, since the energies of different parts of the interface 



add with random signs, i.e. F{x) ^ x^ '■-'^ ^'^^ for x 
in dll), to give 
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M:^L (P-AP, R-AR), 

(3) 



while for F-AF boundary conditions they will add with 
the same sign to give 



(Eint) ^ 



M:$>L (F-AF). 



(4) 



Since the asymptotic forms (||) and (^) are different, it 
follows that the scaling function F{x) in (|l|) will depend 
on the boundary conditions for general x. However, the 
limiting large- a; form, which leads to (^, will be indepen- 
dent of the boundary conditions 

If the limit L M can be achieved in practice, so that 
the form (^ holds, the exponent 9 can be extracted from 
the M-dependence, and is transparently independent of 
the boundary conditions, i.e. if we define 



G(M) = hm L{E,., 

L — >oc 
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then G{M) -> AM'^+'> for M oo {A ^ const.). 
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which is clearly independent of the boundary conditions 
imposed in the x-direction since the limit L — > c» is taken 
before the limit M — > cx). 

In practice we find it convenient to study a broad range 
of R rather than just the regime i? 3> 1, though the ex- 
ponent 6 will ultimately be obtained by extrapolation to 
the R = oo limit. Figure 1 shows a log- log plot of L{Eint) 
against M for various fixed aspect ratios R {1 < R < 32), 
with R-AR boundary conditions. The ground state en- 
ergies were obtained by using exact transfer matrix cal- 
culations, with 2 < M < 12. Each point is an average of 
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FIG. 1. Variation of (Eint) with width M at fixed R for 
R-AR boundary conditions. 



10^ samples. If ARS works perfectly for all L and M, the 
lines corresponding to different R will be parallel, with 
slope 1 + 9. 

The lines are indeed almost straight and parallel. The 
slopes, 9u {R — "random") are presented for different as- 
pect ratios in Table 1. The R ~ 1 result, corresponding 
to squares, is consistent with that obtained earlier using 
the same method However, there is a very slow de- 
crease of the effective exponent with increasing R. We 
have argued that it is sensible to extract the exponent 
from the large- i? limit, since in this limit any errors due 
to finite-size corrections in the a;-direction are eliminated. 
Note that when the data are truly in the regime i? ^ 1, 
the lines in Figure 1 will fall on top of each other. While 
they appear to be approaching a limit, they have not yet 
reached the limit at i? = 32. 

R 6r 6f 

1 -0.289(2) -0.153(2) 

2 -0.286(2) -0.215(2) 
4 -0.285(2) -0.249(2) 
8 -0.283(2) -0.265(2) 
16 -0.285(2) -0.273(2) 
32 -0.283(2) -0.274(3) 

TABLE I. Effective stiffness exponent S as a function of 
the aspect ratio R for R-AR {da) and F-AF (Of) boundary 
conditions, obtained from gradients of the lines in figures hi 
and ^. 

The equivalent results for F-AF boundary conditions 
are presented in Figure 2. In this case the lines are 
clearly not parallel, i.e. there is a strong dependence of 
the measured gradients on R. The lines are also notice- 
ably curved, especially for smaller values of R, e.g. for 
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R = 1 the slope (and therefore the effective value of 9) is 
decreasing with increasing M. 
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FIG. 2. Variation of {Eint) with width M at fixed R for 
F-AF boundary conditions. 



Table 1 contains the effective exponents, Op {F — 
"free"), extracted from Figure 2, for the different values 
of R. These exponents were obtained by fitting straight 
lines to the data, discarding the smallest value of M 
in each case. Note the strong dependence of Op on R. 
For R = \ (i.e. squares) the effective Op differs by al- 
most a factor 2 from the effective Or, but the difference 
gets smaller for larger values of i?, and we argue that it 
approaches zero asymptotically. The argument is based 
on the following observation. For large R, we find that 
the interface energy is nearly always identical, sample by 
sample, for both boundary conditions. The differences 
apparent in Table 1 are due to a small fraction of the sam- 
ples where the energies differ for the two boundary con- 
ditions. In fact, inspection of the interfaces themselves 
shows, as one would expect, that these occupy identical 
locations for the two boundary conditions whenever the 
interface energies are the same. To understand this we 
note the following facts: 

(i) A detailed study of the ground states indicates that, 
for i? ^ 1, the actual spin configurations are independent 
of the boundary conditions in the "interior" regions away 
from the boundaries. The effect of changing the bound- 
ary conditions is localized near the boundary, propagat- 
ing a distance into the system whose mean size is roughly 
proportional to M over the range of M (2 < M < 12) 
studied. 

(ii) The width of the interfacial region also scales roughly 
as M. 

One can understand the latter as follows. Suppose this 
width scales as Af. Then a < 1 would imply that 
the interface is "smooth" on large length scales, in- 



consistent with it having a non-trivial fractal dimen- 
sion ds > d — 1. On the other other hand, a > 1 
would imply that for M ^ 1 the interface is strongly 
"stretched" in the x-dircction. Its energy could then be 
estimated, using (m with M — > M"^ and L — > M as 
Eint ~ M''+('^-i)(^-i)/2 ^ j\,f0^ is inconsistent 

|l2|. We conclude that a — 1. This means that there 
are of order R = L/M independent places in which the 
interface can sit. Any one of these has an energy of order 
, but the prefactor is a random variable with non-zero 
weight at the origin jlSj . The smallest of these therefore 
scales as 1/R, i.e. Eint ~ M^/R for i? ^ 1, a result 
identical to Eq. (|^). From (i) and (ii) we see that the 
probability that the interface enters a region near the 
boundary where the ground-state spin configuration dif- 
fers for the two sets of boundary conditions, and for which 
the interace energies also differ, is of order 1/R for large 
R. 
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FIG. 3. Variation of ^off with system aspect ratio, R, for 
R-AR (upper data) and F-AF (lower data) boundary condi- 
tions. The data are consistent with convergence to a unique 
limiting value, 9 = -0.282(3). 



The above considerations suggest that convergence to 
the limiting behavior (^ occurs as 1/R, so in Figure 3 
we plot the effective exponents listed in Table 1 against 
1/R. The resulting plots are compatible with a linear 
dependence on 1/R, and the data are consistent with 
convergence to a unique value, 6 = —0.282(3), as i? ^ oo. 

It is remarkable that, by exploiting ARS, studies of 
systems with relatively small widths, 2 < M < 12, can 
give results of comparable precision to studies on square 
systems of much larger size. This is especially striking for 
F-AF boundary conditions, where the effective exponent 
depends strongly on the aspect ratio R. For example. Ta- 
ble 1 shows that already for R = 2 the estimate of 9 is 
more accurate than the result 9f = —0.20 obtained for 
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square systems of size up to 24 while for R = 8 
it matches the estimate 9p = —0.266 obtained from 
squares up to size 320 1^. Furthermore, the ARS method 
demonstrates quite convincingly that the stiffness expo- 
nent is boundary-condition independent, a result which 
has not been confirmed conclusively for squares even on 
the largest systems studied This slow convergence 
of 6 for F-AF boundary conditions using squares sug- 
gests that other methods of determining 9, for example 
from the dependence on system size L of the Parisi over- 
lap function P{q) at g = (notably in d = 3), using 
P(0) ~ L^^, as predicted by the droplet model [Q,^, or 
creating droplets by flipping a central spin while holding 
the boundary spins fixed jlj] , might also suffer from large 
finite-size effects. 

We conclude by presenting some results for d — 3. 
The transfer matrix approach restricts us to rather small 
widths, M < 4. For each value of R we obtain an effec- 
tive stiffness exponent O^ff by fitting a straight line to 
the three points M = 2, 3, 4 in a plot of ln(L(|i?[)) against 
InM, and defining the slope to be 1 -I- Oeff- These lines 
are found to be quite straight, so 9eff can be readily ex- 
tracted. It is plotted against in Figure 4, for both 
R — AR and F — AF boundary conditions. Also plotted 
are the equivalent results for the case where free bound- 
aries are used in the directions normal to the domain- 
wall forcing boundary conditions. The final character 
in the legend specifies whether periodic (P) or free (F) 
boundaries have been employed in these directions. The 
free boundary data lie above the corresponding periodic 
boundary data in each case. The dependence of 9eff on 



conditions, the R = 1 result, 9eff — 0.19 is in agreement 
with our earlier results on cubes of side L = 2,3,4 
and those of Hartmann for L < 10 0. For R — > oo, how- 
ever, results for all boundary seem to converge to value 
9 ~ 0.27, significantly larger than previous estimates. 
It is surprising, also, that the boundary conditions in 
the transverse direction do not significantly affect the 
large- i? limit of 9eff- While the small widths M used 
prompts caution in the interpretation of this result (the 
same widths, used for d = 2, would give 9 ~ —0.32 in- 
stead of —0.28), the trends with increasing R are quite 
striking. In particular, the use of F — AF boundary con- 
ditions for cubes (i.e. i? = 1) leads to a very large over- 
estimate of the exponent. 

In summary, the use of aspect-ratio scaling gives, in 
d = 2, results comparable in quality to those obtained 
from square systems of much larger size, and indepen- 
dent of the domain-wall-forcing boundary condition. It 
does this by eliminating finite-size corrections in the di- 
rection normal to the domain wall. For d = 3 the results 
suggest that 9 may be significantly larger than previous 
estimates. 

We thank A. P. Young for a useful discussion. This 
work was stimulated in part by comments from J. M. 
Kosterlitz concerning the possible influence of boundary 
conditions on the determination of stiffness exponents, 
and was supported by EPSRC. 
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FIG. 4. Variation of 9cs with system aspect ratio, R, in 
d = 3 for R-AR (upper data) and F-AF (lower data) bound- 
ary conditions. Periodic and free boundary conditions in the 
y direction are indicated by the final character P or F. 



R in this case is quite striking. For R-AR-P boundary 
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